TAUP 2526-98 
15 October, 1998 



Effect of Non-Orthogonality of Residues in the Wigner-Weisskopf model 

on 

Regeneration for the Neutral K Meson system 
Eli Cohen and L.P. Horwitz* 
Department of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University, Ramat Aviv 69978, Israel 



Abstract: We review the application of the Wigner-Weisskopf model to the two-channel 
decay problem for the neutral K meson system in the resolvent formalism. The residues 
in the pole approximation are not orthogonal, leading to additional interference terms in 
the Ks — Kl, 2tx channel. We show that these terms lead to non-trivial changes in the 
exit beam in comparison to the result calculated with the assumptions of Lee, Oehme and 
Yang, and Wu and Yang, corresponding to semigroup evolution for which the pole residues 
are orthogonal, and hence appear to rule out the applicability of the Wigner-Weisskopf 
model for the computation of regeneration in the neutral K meson system. 
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In 1957 Lee, Oehme and Yang 1 constructed a generalization of the Wigner-Weisskopf 2 
decay model in pole approximation, in terms of a two-by-two non-Hermitian effective 
Hamiltonian, leading to an exact semigroup law of evolution for the two channel K° 
decay. Wu and Yang 3 developed, from this model, an effective parametrization of the 
K° — > 2n decays, resulting in a phenomenology that has been very useful in describing 
the experimental results. Since the early 1970's , Monte Carlo reproduction of the data, 
establishing the parameter values, has been remarkably accurate in the energy 60-120 Gev 
of the kaon beam 3 . 

We wish to point out that the Wigner-Weisskopf model 2 is not, in fact, consistent 
with such semigroup evolution in pole approximation, and that the difference between the 
Wigner-Weisskopf prediction and semigroup evolution could, in principle, be seen experi- 
mentally. The experimental data may therefore rule out the applicability of the Wigner- 
Weisskopf theory even on the level of the pole approximation. In this letter, we consider 
the structure of the two channel Wigner-Weisskopf model on the regeneration process. We 
treat elsewhere the Ks — Kl 2tt interference channel at the output of a regenerator 4 . 

Regeneration systems are designed to reconstitute the Ks beam from a Kl beam in 
the order of K s : K L ^ 10" 3 : 1, since the Kg beam decays rapidly with high branching 
ratio to 27r, while the Kl beam decays to 2tv only on the order of 10 -3 , the measure of 
CP violation. 

The Wigner-Weisskopf model of particle decay (we first treat the one-channel case) 
assumes that there is an initial state in the Hilbert space representation of an unstable 
system; in the course of time, Hamiltonian action evolves this state, and the component 
that remains in the direction is called the survival amplitude: 

A(t) = (^,e- im <f>). (1) 

Defining the reduced resolvent, analytic for z in the upper half plane, in terms of the 
Laplace transform 

POO 

R(z) = -i I exp lzt (0,e- lH V)^ 

(2) 



o 



we see that 



(0, e- lH V) = I dz (</>, -J— 0) , (3) 
Ztvi Jq z — H 



where C is a contour running from oo to slightly above the real axis, and from to 
oo slightly below. The lower contour can be deformed to the negative imaginary axis, 
providing a contribution which is small except near the branch point, and the upper 
contour can be deformed downward (with suitable conditions on the spectral function) 
into the second Riemann sheet, where a resonance pole may appear. Assuming that 
the resonance contribution is large compared to the "background" integrals, the reduced 
motion of the survival amplitude is well-approximated by 

(0,e- iH V) = <7e--°*, (4) 
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where g is the t-independent part of the residue (~ unity), and 

Z = E - 

is the pole position. It appears that this evolution could be generalized for the two channel 
case by replacing zq by an effective 2x2 non-Hermitian Hamiltonian. As we shall see, this 
structure is not generally admitted by the Wigner-Weisskopf model. 

For the two-channel system, consider a state fa, i = 1,2; let us calculate the de- 
cay amplitude fa — ■> |Aj), j = 1,2, the continuum accessible by means of the dynamical 
evolution e~ %m : 



J cU|(A^-^%)| 2 = 1-^1(^,6-* 



Ht \)\ 2 . 



(6) 



Here, the survival amplitude corresponds to the 2x2 matrix 

A ij (t) = (fa,e- iHt fa). 
This amplitude can be approximated in the same way as in (4) by estimating 

(fa,e- im fa) = -L jf ds^-J—^.). ( 7 ) 
It is convenient to write the 2x2 matrix reduced resolvent in the form 

1.x/ 1 



where W(z) is a 2 x 2 matrix. It is almost always true (an exception is the matrix with 
unity in upper right element, and all others zero) that W(z) has the spectral decomposition 

W(z) = gi (z) Wl (z) + g 2 (z)w 2 (z), (9) 

where w\{z) and W2(z) are numerical valued, and gi(z), g 2 (z) are 2x2 matrices with the 
property that 

gi (z)g 2 (z) = 0, (10) 



and 



9\{z) = gi{z), gl(z)=g 2 {z), (11) 



even though W(z) is not Hermitian. These matrices are constructed from the direct 
product of right and left eigenvectors of W(z), and form a complete set 

gi (z)+g 2 (z) = l. (12) 
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Suppose now that the short lived K -decay channel has a pole at zs- The matrix W(zs) 
(in the second Riemann sheet) can be represented as 

W(z s ) = zsgs(zs) + zWs(zs), (12) 

where gsg's = 0, and the corresponding eigenvalues are denoted by zs,z' s . Then the 
reduced resolvent, in the neighborhood of z = zs has the form 

R ij( z ) - 7 vrr, — 77 — ^gs(zs) + t j-rr^- — r ( — ^9s( z s), (13) 

(z - z s )(l + w[{z s )) (z - z' s )(l + w' 2 (z s )) 

where w[(z),w' 2 (z) are the derivatives of the eigenvalues of W(z), of order the square of 
the weak coupling constant (these functions correspond to the mass shifts induced by the 
weak interaction). The residues are therefore gs(zs) an d g's( z s) to a good approximation. 
For the long-lived component, on the other hand, the pole occurs at zl, and at this point, 

W(z L ) = z L9L (z L ) + z' L g' L {z L ). (14) 

The reduced resolvent in this neighborhood is then 

R ij( z ) ~ 7 V7T- I — r ( — ^9l(zl) + t TTTT-T — n — x\9l(?l)- (15) 

(z - z L )(l + w^zl)) (z - z L )(l + w' 2 (z L )) 

Since W(zs) and W(zl) correspond to the matrix-valued function W(z) evaluated at two 
different points, although gs(zs) and g' s (zs) are orthogonal (gs( z s)g's( z s) = 0), in general, 
gs(zs) and ghi^h) are not. If there were no CP violation, the matrix would be diagonal in 
the Ki, K 2 basis, and the two distinct idempotents would be structurally orthogonal even 
though they correspond to the decomposition of the matrix at two different points. One 
can estimate the product 4 ' 5 

upper left element g s g L = 0(a 3 ) = O(10" 4 ), (16) 

in a Lee-Friedrichs type model 6 , where a is the relative amplitude CP violation (indepen- 
dently of the strength of the weak interaction); the other matrix elements are of order a 4 
or a 5 . The regeneration problem was treated previously 5 with emphasis on the question of 
obtaining the maximal possible correction to the semigroup decay law 1 ' 2 . We re-examine 
this problem in the following to establish a lower bound for such corrections, and find that 
the Wigner-Weisskopf predicitons are indeed incompatible with the experimental results. 

We now study the composition of the beam in a regenerator. The usual method for 
the computation of regeneration is not applicable if the U(t) do not form a semi-group. It 
is generally induced by placing a slab of material, such as copper, in the K meson beam. 

At each scattering in the regenerator, the K and K parts of the K$ and K L are 
altered selectively by the amplitudes /, /, the scattering amplitudes for K , K on the 
nuclear target.* After scattering, the new linear combination propagates freely according 
to U'(t). 



* G.Charpak and M.GourdimCERN 67 - 18.0ur / corresponds to 2ni^Atf K \ where 

f K ° is the amplitude used in this reference, and At = ^ Az is the time step of our iterative 
process (Az is the mean lattice spacing of the material). 
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In terms of K and K , the effect of a single scattering is, 

S\K ) = (1 + f)\K ) (17) 
S\K ) = (1+J)\K ) (18) 

where S is the following term: 

S = (1 + /)|K >(^o| + (1 +1)\K )(K \ (19) 
t/ (At) is defined, in pole approximation, in the following way: 

U(At) e" lzsAt (7 S + e- J2iAt (7 L (20) 

where (75, 5^ are the approximate residues of the reduced resolvent, R'(z), on the second 
sheet. We represent these idempotents by 

9s = \K S ){K S \ 
g L = \K L ){K L \ 

where \K$), \K$) are the left and right eigenvectors of the reduced evolution matrix (9) at 
zs for the eigenvalue zs and \Kl), \Kl) at zl for the eigenvalue zl. Let us assume that 
the beam has initially the following form: 

\1>) = a\K ) + p\K ) (21) 

Then, 

U(At)S\iP) = U(At)S(a\K ) + p\K )) = 
= a(l + f)(e- izsAt g s + e-^ At g L )\K ) + 0(1 + J)(e-^ At g s + e~ iz - Ai g L )\K ) (22) 
The following terms appear in this equation: 

<7 S |^o> = \K S )(K S \K ) = * + e ,l = \K S ) 



gs\K ) = \Ks)(K s \K ) = ,\ % = \K S ) 



g L \K ) = \K L )(K L \K ) = J^'t \K L ) 



g L \K ) = \K L )(K L \K ) = " L I = \K L ) 



^(l + lesl 2 ) 

l-e 5 

>/2(l + |cs| 2 ) 
l + e L 



(23) 



Neglecting terms second order in e , / and changing the basis to \K\) , \K2) we get that 
the matrix elements of 

U(At)S = A (24) 



are 

A 11 = (l + f+)e- izsAt 



where 



A12 = (~es + f-)e~ tzsAt + e L e~^ At 
A 2 i = ese~* zsAt + (~e L + f_)e~^ At 
A 22 = (l + f+)e-^ At 

. f + f 

J+ - 



/- = 



2 

f-7 



\K 1 ) = -L(\K ) + \K )) 

\K 2 ) = -L(\K ) - \K )) 

The states \Ki) and \K 2 ) are orthonomal: 

(K 1 \K 1 ) = (K 2 \K 2 ) = 1 



(K 1 \K 2 } = 



We define es, es, e L , e L by 5 : 



\Ks) = -r=±—{\K 1 ) + e s \K 2 )) 

v 1 + ksr 

\Kl)= 7 =±== {e L \K 1 ) + \K 2 )) 

v 1 + M 
fe[= ^_= ((^ 1 | + g s (^ 2 |) 

v 1 + ksr 

fel= _i_= (e L (^i| + (^ 2 |) 



and, 



/ = e i<5 " - 1 



(25) 



(26) 



(27) 



(28) 



/ = e i5 ° - 1 

By multiplying the matrix [U(At)S] n times and keeping only first order in /,/, and e, we 
get the following matrix elements for [U(At)S] n = A(n): 

A(n) u = (l + nf+)e- inzsAt 

A(n) 12 = (i s + f-)e-™ zsAt 

+ e L e~ inZLAt + (e s + + f_)Y n ~ l e-^ kzs+ ^ n - k >^ At 

fc=1 (29) 
A(n) 21 =e s e- tnzsAt + (i L + f-)e- tnZLAt 

+ (es + e L + /_)^™~V l[fczs+(n - fc)2i]A * 
A(n) 22 = (l + nf + )e-^ At 



The product of the residues (16) is 

9s9l = \K s )(K L \(e L + e s ) 
9l9s = \K L )(K s \(e L + e L ) 

and here, if these vanish, 

cl + es = es + £l = 0. 

We obtain for this case, 

A(n)° 12 = (-e L + f-)e~ inzsAt - e s e~ inZLAt + /-V"" 1 e -i[kz s +(n-k)z L ]At 

k=l 

A(n)° 21 = -i Le ~^sAt + (_ eg; + f_y-inz L M + y_y- n - 1 e -i[fc Zs + (n-fc)^]At_ 

We can then substract the two cases, i.e., consider A(n) — A(n)°, and the sums can be 
carried out. We obtain 

( -i(n+2)(z s -z L )At _ -i{z s -z L )At\ -inz L At 

AA(n) 12 = (e L + l s )(e-»— + if >a. _ i ' > 

(30) 

( -i(n+2)(z s -z L )At _ -i{z s -z L )At\ -inz L At V ; 

AA(n hl = (es + iL )(e-<™-' " + ^ >a. _ i ' > 

where A(n)5 2 and ^(n)^ are the matrices corresponding to the orthogonal case. For n 
large compared to unity and we find (the n-dependent exponential factor cancels in the 
ratio), for example, 

AA(n)ia ^ e L + i s 

A(n) 12 (i s + f_) e -*(zs-ZL)At + €l 

The numerator is of order 10~ 4 . Since (zs — z L )At « 1, the denominator is approximately 
eL +es + /-• At IGeV/c, /_ is of order 5 10 -11 , and at higher energies it is smaller (falling 
approximately as /c _1 ); hence ^^^ 2 could be of order unity. The validity of the the 
results of Monte Carlo estimates for regeneration according to the parametrization of ref. 
2 would then appear to rule out the applicability of the Wigner-Weisskopf model for the 
computation of regeneration. 
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